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Abstract 

a rnntroller Performance Evaluation (CPE) 

SSSta. and confirmed dm satisfacior, performance of 
stabilizing ones. 
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stability - 
robustness - 

performance - 


ppfin itioiLS 

all the poles of the system are on the 
left-hand side of the complex plane 
tolerance of system stability to plant 
uncertainty which can be measured m 
terms of minimum singular values or 
gain and phase margins, 
performance of the controller is 
measured in terms of stability and 
robustness. 
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Introduction 

Active controls are becoming an increasingly 
system. 

For examining ihe siabilily «™° n 

ESSSaggs |SH 

touted From these return-difference matrices the 



controller applications including digital flight control 
systems, digitally controlled spacecraft structures, and 
actively controlled wind-tunnel models. 

In the case of actively controlled wind-tunnel models, 
flutter testing adds an extra risk because the controller 
itself can potentially destabilize the model. The on-line 
near-real time CPE methods were employed to check the 
stability of the closed-loop system to reduce the risk of 
damage to the wind-tunnel model and the tunnel. 

The present paper describes the implementation of the 
CPE capability, the structure of the data flow, the signal 
processing methods used to process the data, and the 
software developed to generate the transfer functions. In 
addition, a brief development of the equations used to 
obtain the open-loop plant, controller transfer matrices, 
and return-difference matrices are given. Finally, results 
of applying the CPE methodology to evaluate (in near 
real-time) MIMO digital flutter suppression systems 
being tested on the Rockwell Active Flexible Wing 
(AFW) wind-tunnel model (reference 4) are presented to 
demonstrate the CPE capability. The AFW wind-tunnel 
test program is described in references 5 and 6. 


Controller Performa nce Evaluation 


Closed-loop 

Step 1: Determine the stability margins of the 
closed-loop system during the closed-loop 
testing by evaluating the singular values of 
return-difference matrices, (I+GH) and 
(I+HG). 

Step 2: Determine open-loop plant stability during 
the closed-loop testing to determine the 
open-loop flutter boundary. 


CPE Comnntatinufj 

The CPE computations involve generating frequency 
domain transfer functions of plant outputs, y, due to an 
excitation, u, and control law commands, x, due to the 
same excitation. Fast Fourier Transform (FFT) 
techniques are used to calculate these functions. The 
controller, H, and the return-difference matrices and their 
singular values are then calculated using matrix 
operations. The computations are described in the 
following paragraphs. Figure 2 is a flowchart which 
outlines the CPE procedures. 

Transfer Function* 


A simplified block diagram of the basic closed-loop 
control problem is presented in Figure la. The plant to 
be controlled is represented mathematically by a frequency 
domain transfer matrix, G, with outputs y and inputs e. 
The controller is represented mathematically with a 
transfer matrix, H, with inputs y and outputs x. External 
excitation u is used to excite the system in a specified 
fashion. This excitation is used to derive transfer 
functions between outputs and inputs in either open- or 
closed-loop systems. The open-loop system is one in 
which the control law outputs (commands required for 
controlling plant response) are not fed back into the 
system; i.e., the switch depicted in figure lb is open. 

Controller performance evaluation is a two-mode, 
four-step process. The two modes are open- and closed- 
loop, and each mode can be broken into two steps. The 
process is outlined conceptually for the flutter suppression 
system application as follows: 

Open-loop 


The method used to compute transfer functions is 
described in reference 7. The method was extended in the 
present study to include additional data-windowing 
capabilities and overlap averaging. Windowing 
capabilities include ramp-in/ramp-out, Hanning, cosine 
taper, and cosine bell. The overlap-averaging capability 
allows long time histories to be partitioned into shorter 
time spans, taking advantage of long periods of time 
history data to average out noise. A zero-fill capability is 
available to zero fill time history data to an exact 
increment of a power of two needed for FFT 
computations. The overlap-averaging capability with 
zero-fill provided optimum use of the time history data 
which were obtained. 

The conu-oller-output transfer matrix, X u , is the 
matrix of ratios of the cross-spectra of the controller 
outputs, x, due to the excitations, u, to the auto-spectra of 
the excitations, u; all spectra are obtained from the FFTs 
of the time histories. Each element of the transfer matrix 
is given by: 


Step 1: Verify the controller, H, by comparing with 
the designed control law transfer matrix. 

Step 2: Predict closed-loop performance based on the 
open-loop performance to determine whether 
the control law will stabilize or destabilize 
the system when the loop is closed. 
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matrix, Y u , is given by equauon (2). 

X ( s «iyj (to) ) m 

(y.(«)) u «V: ~~~7~ ' (2) 

X ( Su i u l^)m 

m=l 

The .wo matrices, X„ and Y„. are ihe basis of all .he 
subsequent CPE computations. 

CPF procedures 

rwn i onn Case. Both transfer matrices, X u and 
V 5SS- — system in whichdteloop 

I V - *• — i H£J?K , 52K££ 

r^SS^nalaredennedbyeonanon 

(2), because, in this case, 

G (w) = V u (to). (3) 

The elements of die open-loop 

matrix, HG, (dimensioned n a by n a ) are 

equation (1), since 

HG (cd) = x u («)• ^ 


singular v“SeT!7 KtmSSmahices a. die 
plant input 


Omin (I+HG) 


and the plant output 

o min (I+GH). 

System instabilities occur a. 
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discussing results. 

equations (3) and (5) or (6) since 

GH (©) = G (co) H (co) (7) 

Tbe maximum singular values of the return-difference 
matrices at the plant input, 


<w(i+h<;) 


and (4), 

H x (co) = [Y u (co)Y T u (co )]' 1 Y u (co)X T u (co) < 5 > 

EenXm (6) ^ * 

used. 

H T (co) = [Y T u(co)]" 1 X T u(co) (6) 

To nerform the first step of the open-loop CPE, the 
resulting c ° n “ 0 ^ 

either equauon (5) or ( ) P imp lementation 
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compared for each output/input pair. 

To oerform the second step in the open-loop CPE 

(prcdicitag KnBo! taw will 
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and the plant output 
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system is stable, an encirclement of the critical point 
(the origin) for det(I+HG(<o» indicates that the 
controller is destabilizing. Furthermore, the proximity 
of the determinant locus to the critical point is a direct 
indication of how near the control system is to an 
instability. 

ClQSgd-LQQp CflSC - The objective of step 1 is to 
obtain transfer matrices HG and GH so that singular- 
value, eigenvalue and determinant calculations can be 
performed on the return-difference matrices. The difference 
between closed-loop and open-loop computations is that 
the transfer matrices, X u and Y u , are obtained from the 
closed-loop system. During closed-loop testing, the plant 
transfer matrix is determined from equation (8) 

G T (co) = [\-X^(co)j Y^(<0) (8) 

and the open-loop controller-plant transfer matrix is 
determined from equation (9) 

[HG] T (to) = J^I-X^((o)j xj(«). (9) 

Noting that the quantity 

['-*>)] ' 

is in both equations, transfer matrices G and HG are 
obtained simultaneously using matrix partitioning. The 
controller transfer matrix, H is calculated using either 
equation (5) or equation (6), as in the open-loop case. 

The singular-value plots and eigenvalue plots are 
interpreted the same way for closed-loop testing as they 
were for open-loop testing. However, care must be taken 
in interpreting the determinant plots. If the open-loop 
plant is unstable with one right-hand pole (for positive 
frequencies), then the plot will show no net encirclement 
about the origin for a stable closed-loop system. As with 
the open-loop case, the proximity of the determinant 
locus to the critical point and the proximity of the 
minimum singular values to zero are used as measures of 
closed-loop stability. 

To obtain results for step 2 of the closed-loop mode, 
the inverted maximum singular value of the plant transfer 
matrix, G, as given by, 

1 

<W{G(co)) 

was an excellent indicator of poles in the proximity of the 
imaginary axis. The frequency of instability was 
determined by where the inverse maximum singular value 
approaches zero. Tracing values of closest approach was a 


useful way of determining the open-loop plant flutter 
boundary with respect to some changing test condition, 
such as dynamic pressure. 


Summary ftf Fllltter-Sunnression Testing 

During flutter suppression testing, the control 
systems were operated in both open-loop and closed-loop 
modes. For the purpose of maintaining both model and 
tunnel safety, each candidate control law was initially 
tested open loop to insure that the control law itself 
would not destabilize the wind-tunnel model during 
closed-loop tests. The feedback was digitally switched 
open at the control law output point and the response data 
was collected at the control law input and the control law 
output locations. The stored data are the excitation, u the 
control law output, x, and the control law input, y, 
depicted in figure 1. The appropriate transfer functions 
were generated from these responses and then the CPE 
capability was exercised to predict the stability of the 
closed-loop system (figure 2). If the control law was 
predicted to be stable, the switch was closed and the 
closed-loop flutter suppression (FS) testing for that 
candidate control law commenced. During the closed-loop 
portion of the test, the excitations were inserted and 
response data was taken at the same test conditions as the 
initial open-loop tests. At each test point, stability 
margins and open-loop plant stability were determined 
before proceeding to the next test point. 

Figure 4 shows the plant and FS control laws. For 
this model, separate FS control laws were required for 
controlling symmetric and antisymmetric flutter modes 
simultaneously. The CPE excitation was input to the 
control surfaces either symmetrically or antisymmetrically 
depending on which control law was being evaluated. The 
responses, y, were summed or differenced before saving. 


Inscription nf the CPE implement^,, 
Hardware 

The digital controller not only provided basic control 
and flutter suppression of the model, but was the source 
for the excitations needed for CPE. The excitation was 
generated digitally and added to the control law actuator 
commands. The digital excitation, actuator commands, 
and sensor measurements used by the control law were 
stored and then transferred to a SUN 3/160 computer with 
a SKY Warrior II array processor board (capable of up to 
8K FFT compulations) where the FFT computations were 
performed. The matrix computations to obtain the 
singular values of return-difference matrices were also 
performed on the SUN 3/160. 

Software 

All the digital controller software was written in the 
high level C programming language, except for the 
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commands required to perform the actual Rations on 
the array processor. Operation code command blocks were 
generated for these. The FFTs of the time histones o 
the excitations and the responses, and the transfer 
tacSs computed in a Fowun-77 program 

optimized to take advantage of the ^ t ° r ' pr ? c ® s *p| 
capabilities on the array processor. The detailed CPE 
analysis capability was implemented using MATLAB 
software operations (reference 8). 

PiHillltri and discussion 

Both SISO and MIMO flutter suppression control 
laws were designed for the AFW wind-tunnel model 
During the wind-tunnel test the FS control laws wer 
successfully tested and their performance evaluated using 
the CPE capability presented in this paper. 

The data for performing CPE was obtained by 
exciting, one at a time, all pairs of control surfaces used 
by the control law. The surfaces were excited 
symmetrically or antisymmetrically (correspond.ngto 
either the symmetric or the antisymmetric FS contro 
laws) with a logarithmic sine sweep. In most cases the 
sJeep duration was 150 seconds over a frequency range 
from 4 to 35 Hz. The excitation, the plant responses as 
measured by various accelerometers, and the controller 
outputs, corresponding to actuator commands, weresaved. 
While these stability computations were being performed 
for the symmetric controller, the antisymmetric 
excitations were performed and the transfer matrices 
generated. Final results and plots were available 
approximately 3 minutes after the last sweep was 
performed. 

The first five seconds of a sample excitation, 
resulting actuator commands and accelerometer responses 
are shown in figure 5. The quantization levels of the 
digital controller can be noted by the step-like behavior of 
the time histories for the actuator commands. A typical 
transfer function at one test condition is shown in figure 

6 . 
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obtained from the universal gain and phase diagram of 
figure 3. Since the a min = 0.37, the gain m^gmforzero 
phase margins are approximately -2.7 db and +4 0 db 
for a 20 degree phase margin, the gain margins are -0. 

and +2.05 db. 

MIMO Co ntrol Law 

Results for a MIMO FS control law design 
(control law 1 of reference 6) are presented next. Since 
the plots of the maximum and minimum singular values 
of the return-difference matrices at the ^ J 

output locations along with the c.genvalueplot fo Ui 
MIMO system are not identical as in the SISO 
system, all three plots are shown in figure 9. In the 
upper right plot of figure 9, the ratio of the maximum to 

theminimum singular value (ie c ° ndltlon /^ s J 
large in the vicinity of 10 Hz. As discussed previously 
large condition numbers indicate uncertainty in the 
computation of the minimum singular valueS m J he 
determinant plot (lower right) shows an encirclement of 
the origin where the open-loop plant is known to be 
stable, thus indicating that the MIMO control law wo 
be destabilizing. Both the minimum singular value and 
* g Sue plot show noticeably low stab.li.y margin 
near 20 Hz. Upon further investigation, one of the 
elements of this desUibilizing controller transfer matrix 
was examined and is shown in figure S 

shows a peak magnitude close to 20 Hz. 

FlUttff Prediction 

A plot of the inverse maximum singular values 
of the plant transfer matrix computed from closed-loop 
tests for a dynamic pressure of 200 psf is shown in figure 
11 The frequencies at which the two local minima 
(approximately 8.5 and 11 Hz) of >nverse max,mum 
singular value curve occur are indicated by the arrows 
These frequencies correspond to the frequencies of th 
modes which coalesce to create flutter. These frequencies 
were determined for many dynamic pressures and figure 12 
shows a plot of the two frequencies as a function of 
dynamic pressure. The dynamic pressure at whic^te two 
curves appear to coalesce indicates a potential point for 
open-loop flutter. 


Typical CPE results for a symmetric SISO control 
law (control law 2 of reference 6) obtained during the 
closed-loop wind-tunnel tests are shown in figures 7 and 
8 Because, for SISO control laws, eigenvalue and 
singular value plots are identical, only one plot is shown. 
Also for SISO control laws, the determinant plot is 
identical to the Nyquist plot. The i determinant plot in 
figure 7 shows no encirclement about the origin i (Ore 
critical point) at a dynamic pressure of 200 psf where the 
open-loop plant is known to be stable. Figure 8 show 
the CPE results of a closed-loop system where the plant 
is unstable. Since the minimum singular values _are 
greater than zero and there is no net encirclement of the 
critical point, these plots indicate that the controller is 
‘SS plant. Using the minimum stngnla, value 
from figure 7, guaranteed stability margins can be 


P onelnsionS 

A Controller Performance Evaluation (CPE) 
methodology was developed to evaluate the performance 
of multivariable, digital control systems. The method 
was used and subsequently validated during the 
tunnel testing of an aeroclastic model equipped with a 
digital flutter suppression controller. Through the 
effort a wide range of sophisticated real-time analysis 
tools were developed. These tools proved extremely useful 
and worked very well during wind-tunnel l ^ sUn S- 
Moreover, results from open-loop CPE were t e so e 
criteria for beginning closed-loop testing. I" this way, 
CPE identified potentially destabilizing controllers before 
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actually closing the loop on the control system, thereby 
avoiding catastrophic damage to the wind-tunnel model or 
the tunnel. Open-loop plant transfer functions derived 
from CPE computations were used to redesign and 
improve control laws. CPE results also proved useful in 
determining open-loop plant stability during closed-loop 
test conditions. F 
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a) Basic closed-loop control problem. 



b) Open-loop control problem. 

Figure 1. Controller-Plant diagrams depicting the 
control problem. 
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Figure 2. Flowchart of CPE procedures. 




Figure 4. Flutter control law implementation. 
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Figure 5. Time histories of sample excitation, actuator 
commands and accelerometer responses. 


7 










10 


a(l+HG) 


0 

u Freq(Hz) 40 


a mjn =0.43 at 6.4 Hz 



a) Singular value plot. 
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Figure 8. Closed-loop CPE results for a symmetric 
SISO control law (open-loop plant is 
unstable), M=0.44, q=250 psf. 


Figure 7. Closed-loop CPE results for a symmetric 

SISO control law (open-loop plant is stable), 
M=0.39, q=200 psf. 
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Figure 1 1 . Plot of inverse maximum singular values of 
the symmetric plant transfer matrix, q 
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Figure 12. Frequencies of modes contributing to 

symmetric flutter determined from inverse 
of maximum singular values. 
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